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Abstract 



The effects of radial advection and thermal diffusion were considered in investigating the 
linear stability of an optically thin, two-temperature accretion disc. If the disc has only 
very little advection, we proved that the thermal instability exists when the disc is geo- 
metrically thin. But it dispears in a geometrically slim disc if the thermal diffusion was 
considered. Moreover, if the disc is advection dominated, the thermal instability does not 
exist. In addition, we found that the instabilities of inertial- acoustic modes exist only in 
a geometrically thin disc or an advection-dominated disc with low Mach number, whereas 
the Lightman & Eardley viscous instability always dispears in a two-temperature disc. 
A simple comparison also showed that an optically thin, bremsstrahlung cooling domi- 
nated disc is generally more thermally unstable than a two-temperature disc if it is not 
advection-dominated. 

Keywords: accretion, accretion discs - black hole physics - instabilities 
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1 Introduction 



The optically thin, two-temperature accretion disc model was suggested in 1970s to ex- 
plain the hard X-ray spectra observed in black hole candidates such as Cyg X-l (Shapiro, 
Lightman & Eardley 1976, hereafter SLE). The standard geometrically thin, optically thick 
accretion disc model, proposed by Shakura & Sunyaev (1973), was unable to account for 
them. In a two-temperature disc, the temperature of electron is about 10 9 K and the ion 
temperature is one or two orders higher. Therefore, the two-temperature disc has much 
higher temperature than the cool, optically thick disc and can produce the observed hard 
X-ray spectrum above ~8 keV. In last two decades, SLE model has been widely stud- 
ied and applied in modeling of X-ray binaries and active galactic nuclei(e.g., Kusunose 
& Takahara 1988, 1989; White & Lightman 1989; Wandel & Liang 1991; Luo & Liang 
1994). However, the early work of Pringle (1976) and Piran (1978) indicated that the two- 
temperature disc model is still thermally unstable. As in the standard thin disc model, 
which is both thermally and viscously unstable in the inner region (Shakura & Sunyaev 
1976; Lightman & Eardley 1974), the rapid growing of instability may result in the break 
down of the disc equilibrium and make it unlikely to be the real configuration of accretion 
flow. 

In most previous studies we mentioned above, the disc model was relatively simple and 
many effects were ignored. Recently, the effects of radial advection were extensively studied 
in both optically thick and optically thin accretion discs (Abramowicz et al. 1988; Kato, 
Honma & Matsumoto 1988; Narayan & Yi 1994, 1995a,b; Abramowicz et al. 1995; Chen et 
al. 1995; Chen 1995; Misra & Melia 1996; Nakamura et al. 1996). If the radiative cooling 
is not efficient to balance the viscous heating in the disc, some energy will be advected 
inward and the advection will be not negligible. It has been suggested both the disc 
structure and the stability properties of an accretion disc with advection term included 
are different from those in previous models where advection was totally neglected. By 
analyzing the M(E) slope of disc structure and comparing the cooling and heating rates 
near each equilibrium curve, Chen et al. (1995), Abramowicz et al. (1995) and Narayan 
& Yi (1995b) have suggested that an advection-dominated disc is both thermally and 
viscously stable whether the disc is optically thin or optically thick. However, such a 
stability analysis, as well as Piran's criteria, can only be applied to the case with long- 
wavelength perturbations (Chen 1996). More recently, Kato, Abramowicz & Chen (1996) 
performed an analytic stability analysis of the advection-dominated discs by considering 
the short wavelength perturbations. They indicated that the advection-dominated disc is 
thermally stable if it is optically thin but thermally unstable if it is optically thick. This 
result was confirmed by a subsequent study of Wu & Li (1996), who considered not only 
the stability of thermal mode, but also the stability of viscous mode and inertial-acoustic 
modes. They also concluded that the thermal diffusion has a significant contribution to 
stabilize an optically thin, advection-dominated disc but enhance the thermal instability 
of an optically thick, advection-dominated disc. 

Although some authors have mentioned that an optically thin, two-temperature disc 
may be thermally stable if it is advection-dominated (e.g., Narayan & Yi 1995b), however, 
no detailed stability analysis has been done so far to consider the effects of advection on the 
stability of a two-temperature accretion disc. For example, in the recent work of Wu & Li 
(1996), only the optically thin disc with bremsstrahlung radiative cooling was addressed. 
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We noted that advection may be quite important in a two-temperature disc as calculated 
by some authors (for example, see Esin et al. 1996). Moreover, a hot two-temperature 
disc is probably not geometrically thin and the Piran's criteria may not be applicable in 
this case. Therefore, in this paper we performed a detailed linear study to investigate 
the stability of a hot two-temperature disc, considering both the effects of advection and 
thermal diffusion. Besides the stability of thermal mode, the stability of viscous mode and 
that of inertial-acoustic modes have also been investigated. 



2 Basic equations 

The simplified equations which describe the equilibrium structure of a two-temperature 
disc has been given by SLE. However, in order to study the stability properties, we must 
involve the detailed time- dependent hydrodynamic equations. In addition, because we 
will investigate the influences of some terms such as the radial viscous force, advection 
and thermal diffusion on the disc stability, we should also included them in the equations. 
Ignoring the self-gravitation and adopting the pseudo-Newtonian potential (Paczyhski & 
Wiita 1980), ^ = —GM/(R — r g ), we can write the vertical integrated time- dependent 
equations in a cylindrical system of coordinates as follows: 

s f + EV; f -=(n'-tUV = -2^ + F„ (2) 

&3 an + SrV , djrHti = 3_ { ^ r3 9n (3) 

dt r dr dr dr ' 

c » |E f + Ev; f " (rs " 1)r( f + v; f >i = E " (r f » 2 " Q - + Q " (4) 

where V r , Q are the radial velocity and angular velocity, Qk is the Keplerian angular 
velocity, given by ilk 2 = (f§^)z=o- P, T, £ and H are the total pressure, temperature, 
surface density and the disc height, C v and r% are the heat capacity per unit mass and a 
quantity associated with f3, the ratio of gas to total pressure. For an optically thin, two- 
temperature disc, (5 — 1, C v — 3p/ (2pT), T 3 = 5/3 and p = pk{T^ + T e )/m p , where T iy T e 
and m p are the ion, electron temperatures and proton mass. F u is the radial viscous force, 
which is often neglected in geometrically thin accretion discs but is perhaps not negligible 
and should be considered in accretion discs with advection (Narayan & Yi 1995a). It is 
given by (Papaloizou & Stanley 1986) 

u ~ dr [ 2 r dr 1 r dr ' W 

where v T is the kinematic viscosity acting in the radial direction. In this paper we take 
v r = u, where v is the viscosity acting in the azimuthal direction and is expressed as the 
standard a prescription (Shakura & Sunyaev 1973), v = ac s H. c s is the local sound speed 
defined by c s = p/p. The disc height is given by H — c s /Qk- Q- at the right side of 
Eq. (4) represents the radiative cooling. For an optically thin, one-temperature disc, the 
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radiative cooling mechanism is usually taken as thermal bremsstrahlung. However, for a 
two-temperature disc, the cooling mechanism is unsaturated Comptonization through the 
loss of energy of electrons. In this case, the cooling rate can be expressed as (see also SLE) 



Q_ = (4kT e / m e c 2 ) pH K es ll r c 



(6) 



where m e is the electron mass, K es = 0Acm 2 g~ 1 is the electron scattering opacity. U r is the 
radiation energy density of soft photons, which we assume, for simplicity, does not change 
on timescale short compared to (Pringle 1976). Because the ions and electrons are 
coupled by collisional energy exchange, the loss of energy of electrons can be balanced by 
the energy capture from ions. The exchange rate is 



where ve is the electron-ion coupling rate, and can be approximated by ve = 2.4 x 
10 21 ZnApT e ~ 3 / 2 (Spitzer 1962) where the Coulomb logarithm InA is about 15 (SLE). Taking 
Q_ = 4- and % » T e , we can get Q„ oc S 7 / 5 T/ /5 . Q t at the right side of Eq. (4) 
represents the thermal diffusion defined as Q t = V ■ (KVT), where K is the vertical 
integrated thermal conductivity given by K = HC v v = afil k H 3 p/T, where / = 3(8—7/3)/* 
and /* is of the order of unity (Kato et al. 1996; Wu & Li 1996). In this paper, we take 
/ = 3 since a two-temperature disc is usually gas pressure dominated and /3 ~ 1. 

The equilibria of a two-temperature disc can be calculated by numerically solving 
above equations with d/dt=0. Ignoring the influences of radial viscous force and thermal 
diffusion, the disc equilibria have been calculated by many authors (e.g., Narayan & Yi 
1995b; Chen et al. 1995; Nakamura et al. 1996). In this paper, we assume the influences 
of radial viscous force and thermal diffusion on the disc equilibria are small, and refer to 
the solution of the disc structure obtained by Narayan & Yi (1995b). We also restrict 
our stability analyses within the validity of local approximation and vertical integrated 
equations, which requires that A < r and kV r < Qk- These two inequalities can be 
summerized as jj > A > 2na^, which is well satisfied for a geometrically thin accretion 
disc in a wide range of perturbation wavelength but for a geometrically slim disc only when 
the perturbation wavelength is short and viscosity efficient a is sufficiently small (Kato et 
al., 1996; Wu & Li 1996). The vertical hydrostatic balance can be realized if the time scale 
associated with the perturbations is longer than the dynamical time scale. Although for a 
geometrically thick disc, the perturbation wavelength can be smaller than the disc height 
(A < H), which means that the perturbations are also local in the vertical direction (Kato 
et al. 1997), the consideration of the local vertical perturbations for a two temperature 
disc does not change the stability properties very much (Yamasaki 1997). 

3 Linear perturbations and dispersion relation 

In this section, we consider the linear perturbations to the two-temperature disc. The 
radial perturbations of V r , f2, £ and T are assumed of the form (SV r , Sfl, ST) ~ e »(wt-fcr) ^ 
where k is the perturbation wavenumber defined by k = 2tt/X, A is the perturbation 
wavelength. Taking the perturbed quantities to the basic equations given in last section 
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(3/2)pHv E k(T l - T e )/m. 



(7) 
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and considering the local approximation, we can obtain following perturbed equations: 

_ 5T. . e 5V r „ 

^^-^7^=0, 8 
S H9 k r ' 1 ; 

~<5E . 4 2N 5\/ r -50 ~<ST 



xaeHg— + 4,—- + (a + ae 2 ) — + laeHg— = 0, (10) 
E 20 fc r 9 k T 



2 7 2 5H ag# 2 5K r 2iae# £0 ,,3 ,4 1 2 2 ^ T n 
" ((7 + a( -5 + 5^ } ^-^F?V + ^-?4 +( 2 a + a( 5 + 5^ } Y" = °' (U) 

where a = cr/9 k , an d cr = i(u — kV r ). 9 = 9/9 k , H = H/r, and e = kH. g = 4 — 29, 
and x — x/^k where x is the epicyclic frequency defined by x 2 — 20(20 + r|p). m is 
the Mach number defined by m =\ V r \ jc s . q is the ratio of advective energy to viscous 
dissipated energy, namely 

aPK§ - (T, - l)T(§ + K§ )] = <,SK-f) 2 . (12) 

If the disc is radiative cooling dominated, q is nearly zero and if it is advection-dominated, 
q is nearly 1. 

By setting the determinants of the coefficients in above perturbed equations to zero, 
we get a dispersion relation: 

aid- 4 + a 2 cx 3 + a 3 <7 2 + a 4 a + a 5 = 0, (13) 

where a,(i = 1, 5) is the coefficients given by 

3 

ai = 2> 

«2 = «[e 2 (/ + ^)-(^ + ^ 2 ], 

a 3 = aefl( 2 U e (^ Z_o) _ i-ll + I( ae 2 ) 2 (6 + 7/) + -e 2 + -y 2 , 

y L V 15 15 y; m J 3 V ; V JJ 2 2 A ' 

a 4 = 2ia 2 9eg 3 - + ag 2 [-U\ + ^g)(ae 2 ) 2 - iae 3 ^- + ^(ae 2 ) 2 - (- + \q)x 2 
m 355 m 3 55 

+j|(<? - l)e 2 ] + <x?g[t - 50] + ae 2 [^/(« e 2 ) 2 + fx 2 + ^ + /)], 

a 5 = («e) 2 [y (1 - qW + *V(-f + p) + ^("| - \<1 ~ 2ffl) + e 4 /]- 

The stability properties of two inertial-acoustic modes, thermal and viscous modes can be 
obtained by analyzing the four kinds of solutions of the dispersion relation. The real parts 
of these solutions correspond to the growth rates of the perturbation modes and the imag- 
inary parts correspond to their propagating properties. If the growth rate corresponding 
to certain mode is positive, this mode is unstable. Otherwise, it will be stable. Among 
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the four kinds of modes in the disc, the stability of thermal and viscous modes have been 
more extensively studied in the literatures. The instability of inertial-acoustic mode (or 
called the pulsational instability), was first addressed in viscous accretion discs by Kato 
(1978) and was detailed studied by Blumenthal, Yang & Lin (1984) later. It has been sug- 
gested that the inertial-acoustic instabilities may account for the observed quasi-period 
oscillations (QPO) in some Galactic black hole candidates (Chen & Taam 1995). 

4 Numerical results 

In this section we will numerically solve the dispersion relation Eq. (13). The rotation 
law of disc is further assumed to be Keplerian. Although Narayan & Yi (1994) obtained 
a sub-Keplerian rotation law (fl ~ 0.5Qk in a hot optically thin disc), such a difference 
will not affect significantly on our results. According to the different possible structure of 
a two-temperature disc, we solve the dispersion relation in following three cases: 
(a) Geometrically thin, cooling dominated two-temperature disc. We take f2 = x = 1, 
a = 0.01, m = 0.01. According to the local restrictions, X/H is set from 1 to 80 for a 
geometrically thin disc (in this section we take H/r = 0.01). By solving the dispersion 
relation, we get the results shown in Fig.l. For comparison, we also show the stability 
properties of an optically thin, geometrically thin bremsstrahlung disc, which has been 
suggested to be thermally unstable (Pringle, Rees & Pacholczyk 1973). Fig. 1(a) shows 
the case without advection. We can clearly see that the thermal mode is always unsta- 
ble and the Lightman & Eardley viscous mode is always stable. This result agree well 
with the previous finding. The inertial-acoustic modes are slightly unstable but stable 
to very short wavelength perturbations. In addition, we see that the thermal mode in 
a bremsstrahlung disc is more unstable than that in a two-temperature disc. Fig. 1(b) 
shows the case with very little advection (q = 0.01). It is clear that the inclusion of very 
little advection has nearly no effects on the thermal and viscous modes, but leads to the 
departure of two acoustic modes. Comparing to the case without advection, the outward 
propagating acoustic mode (hereafter O-mode) now becomes more unstable, whereas the 
inward propagating acoustic mode (hereafter I-mode) becomes more stable. In Fig. 1(c) 
we show the case with q = 0.01 and m = 0.1. Comparing with Fig. 1(b), the increase 
of Mach number leads to the decrease of the departure of two acoustic modes. Actually, 
as we have noted in a previous work (Wu & Li 1996), such a departure is proportional to 
the term q/m. Both the O-mode and the I-mode are unstable to the longer wavelength 
perturbation if q/m less than 0.1. The change of q/m, however, has nearly no effects on 
the thermal and viscous modes. Fig. 1(d) shows the case when q = 0.01, m = 0.01 but 
with thermal diffusion included. Comparing with Fig l.(b), we see that the inclusion of 
thermal diffusion has nearly no effect on the viscous mode and two acoustic modes, but 
it tends to stabilize the thermal instability in the short perturbation wavelength case. In 
general, if the disc is geometrically thin and cooling dominated, the stability properties of 
a two-temperature disc are quite similar to those of a bremsstrahlung disc. But we still 
see clearly that the thermal mode in a bremsstrahlung disc is more unstable than that in 
a two-temperature disc. We noted that these results agree well with those obtained by 
Luo & Liang (1994), who have pointed out that the thermal mode is always unstable and 
viscous modes is always stable in a hot optically thin disc. However, they did not consider 
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the stability of inertial-acoustic modes and the effects of radial viscous force, advection 
and thermal diffusion. 

(b) Geometrically slim, cooling dominated two-temperature disc. Some previous works have 
suggested that a hot optically thin disc may be not geometrically thin but geometrically 
slim or thick (e.g., SLE) even if the advection term was ignored. Here we show the stability 
of a geometrically slim, two-temperature disc. In this section we take H/r = 0.6. The 
perturbation wavelength is set from 0.2 to 2 according to the local approximation. Fig. 
2(a) shows the case when m = 0.01, a = 0.001 and without advection (q = 0). We can 
see that the inertial-acoustic modes and the viscous mode are always stable. Only the 
thermal mode is unstable when perturbation wavelength is longer than 1.5H. Fig. 2(b) 
show the similar results as Fig. 2(a) but with very little advection included (q = 0.01). 
In comparison with Fig. 2(a) we see that two acoustic modes now slightly depart from 
each other. For clearity, Fig. 2(c) compares the stability of thermal modes in two cases 
above. The stability of thermal mode of an optically thin, bremsstrahlung disc is also 
shown for comparison. It is quite clear that the inclusion of very little advection has only 
very slight effects on the thermal mode, which is always unstable in a bremsstrahlung disc 
but can become stable in a two-temperature disc if the perturbation wavelength is shorter. 
The case with thermal diffusion considered is shown in Fig. 2(d). We clearly see that all 
four kinds of modes in a two-temperature disc are always stable if the thermal diffusion 
is considered. We also noted that an optically thin, geometrically slim bremsstrahlung 
disc can become stable if the thermal diffusion is considered, even if it is more thermally 
unstable than a two-temperature disc when the thermal diffusion is ignored. 

(c) Geometrically slim, advection- dominated two temperature disc. The advection-dominated 
equilibrium of an optically thin, two-temperature disc has been recently constructed by 
some authors (e.g., Narayan & Yi 1995b; Chen et al. 1995; Nakamura et al. 1996). Such a 
disc is usually geometrically slim. In Fig. 3, we show its stability against short wavelength 
perturbations. Here we take q = 0.99 and H/r = 0.6. Fig. 3(a) shows the case with 
m = 0.01 and a = 0.001. In this case, the thermal, viscous modes are always stable. 
The inertial-acoustic modes depart from each other. The I-mode is stable but the O-mode 
can become unstable if the perturbation wavelength is larger. Fig. 3(b) shows the case 
with different Mach number (m = 0.1). Comparing with Fig. 3(a), the departure of two 
acoustic modes becomes less due to the decrease of q/m. The O-mode can become stable 
if q/m less than 10. However, the change of q/m has no effects on the thermal and viscous 
modes, which are always stable. In Fig. 3(c), we show the case similar to Fig. 3(a) but 
with thermal diffusion included. We can clearly see that the inclusion of thermal diffusion 
has a significant effect to stabilize the thermal mode, but it has nearly no effects on the 
viscous and acoustic modes. We noted that in an optically thin, advection-dominated disc, 
the stability properties do not change evidently if the different radiative cooling mecha- 
nism was involved. Whether the disc is bremsstrahlung one or two-temperature one, the 
stability properties are always the same. This is, of course, due to the unimportance of 
radiative cooling in an advection-dominated disc. 
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5 Discussions 



The stability of a hot optically thin two-temperature disc has been discussed in above sec- 
tions by numerically solving the dispersion relation. We proved that the thermal instability 
exists in the disc if it is cooling dominated but dispears if it is advection-dominated. The 
Light man- Ear dley viscous instability is always absent in a hot optically thin disc. These 
properties agree well with previous qualitative results. We have also investigated the sta- 
bility of two inertial-acoustic modes, which has not been discussed previously in a hot 
optically thin, two-temperature disc. In a geometrically thin, radiative cooling dominated 
disc, we found that the O-mode is always unstable but the I-mode is unstable only if the 
term q/m less than 0.1. If the disc is cooling dominated but geometrically slim, no inertial- 
acoustic mode is unstable. However, in a geometrically slim and advection-dominated disc, 
the O-mode can become unstable if q/m larger than 10, though the I-mode is always sta- 
ble. These stability properties of acoustic modes may be important when they are involved 
to explain the QPO phenomena in some systems such as Galactic black hole candidates 
(Chen & Taam 1995; Manmoto et al. 1996). 

A simple comparison shows that the bremsstrahlung disc is more thermally unsta- 
ble than the two-temperature disc if they are cooling dominated. This can been seen 
clearly by comparing their different temperature-dependence of the cooling rates. In a 
bremsstrahlung disc, the cooling rate Qb~ is proportional to S 2 . In a two-temperature 
disc, the cooling rate Q t ~ is proportional to E^T 1 / 5 . In these two cases the energy gen- 
erated rate by viscous dissipation can be both expressed as Q + = Xz/f2|^ oc ST. Thus, 
if we give a positive perturbation 5T and keep the S unperturbed, Q + will grow more 
rapidly than Q b ~ and Q t ~ , which leads to thermal instability (see also Piran 1978). Be- 
cause Q + — Q t ~ is less than Q + — Qb~ with a positive perturbation of temperature, the 
two-temperature disc is less thermally unstable than the bremsstrahlung disc. If the hot 
optically thin disc is advection-dominated, however, the disc stability will not depend on 
the detailed cooling mechanism due to its less importance. 

It is quite evident that the advection and thermal diffusion have significant effects on 
the stability of a hot optically thin disc. These effects are more evident if the disc is not 
geometrically thin. As having been pointed out by Wu & Li (1996), the thermal diffusion 
in an accretion disc is in proportional to (H/r) 2 . Therefore, the thermal diffusion should 
be seriously considered in the hot optically thin disc which is usually geometrically slim 
or thick. We think that the inclusion of thermal diffusion will affect not only the stability 
but also the structure of accretion discs. A future investigation on this point is expected. 

More recently, we noted that Kato et al. (1997) pointed out that the perturbations 
could be also local in the vertical direction for a geometrically thick disc. This is true but 
the consideration of the vertical local perturbations seems to have only minor contribution 
to the growth rate of unstable mode. For a two temperature disc, the independent analytic 
stability analysis of Yamasaki (1997) found that there are two kinds of modes in the disc. 
One is thermal mode which is slightly unstable and the other is viscous Lightman & Eardley 
mode which is always stable. The inclusion of thermal diffusion will stabilize the thermal 
mode. These result agree quite well as our quantitative analyses. However, we noted that 
the radial viscous force was not included in the analysis of Yamasaki (1997). Consireation 
of this will contribute a term in proportional to (kH) 2 in the perturbed equation of radial 
momentum conservation, which can not be neglected for a geometrically thick disc. It will 
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result in the absence of the weak thermal instability of a hot optically thin disc when the 
short wavelength perturbations are considered (Wu, Yang & Yang 1994; Wu & Li 1996). 

We should noted that the radiative cooling rate of a two-temperature disc may be not 
as simple as we assumed in Section 2. For example, we have assumed that the radiation 
energy density of soft photons, U r , is unchanged on dynamical timescale, which may be 
not always the case. Some recent works on the hot optically thin discs suggested that the 
cooling mechanism is rather complicated (See e.g., Esin et al. 1996). The disc structure 
also depends on the radius significantly. This implies that the disc stability may be quite 
different from one radius to another. Thus, the global stability analyses of a hot optically 
thin disc are still needed. In addition, the pair production and annihilation in a hot 
accretion disc are totally neglected in this paper for simplicity because the inclusion of 
pair process will make the linear stability analysis very complicated. Although this process 
may be important especially if the disc temperature is higher enough, the pair density is 
believed to be not too high and their influence is limited (Bjornsson et al. 1996; Kusunose 
& Mineshige 1996). However, we still expect a future detailed work would be done to 
investigate their influences on the structure and stability of a hot accretion disc. 

Together with some previous detailed stability analyses, our study show that there 
are perhaps only two stable thermal equilibria of accretion discs. One is optically thin, 
advection-dominated and the other is optically thick, gas-pressure and radiative cooling 
dominated. These equilibria are probably related with some stable, inactive astronomical 
systems. For example, the optically thin, advection-dominated accretion discs have been 
suggested to exist at the center of our Galaxy (Narayan, Yi & Mahadevan 1995), nearby 
elliptical galaxies (Fabian & Rees 1995), fainter AGNs (Lasota et al. 1996) and soft X- 
ray transient sources in the quiescent state (Narayan, McClintock & Yi 1996). On the 
other hand, the unstable thermal equilibria of accretion discs, such as the optically thin, 
radiative cooling dominated one, the optically thick, radiation pressure dominated one 
and the optically thick, advection-dominated one, probably exist in some unstable systems 
such as the inner region of AGNs, X-ray binaries and cataclysmic variables. In order to 
understand the light variability of these systems more clearly, the detailed time-dependent 
nonlinear studies on the evolution of their disc structures are still expected. 
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FIGURE CAPTIONS 



Figure 1. Stability of a geometrically thin, optically thin disc. Heavy lines represent the 
case of a two-temperature disc and light lines represent the case of bremsstrahlung disc. 
In (a) no advection is included, the solid, long-dash and short-dash lines correspond to 
the acoustic modes, thermal mode and viscous mode respectively, (b), (c) and (d) show 
the cases with very little advection (q=0.01), parameters (m, f) are taken as (0.01, 0) in 
(b), (0.1, 0) in (c) and (0.01, 3) in (d). The long-dashed and short-dashed lines in (b), (c) 
and (d) have the same meaning as in (a) but the solid and dotted lines correspond to the 
O-mode and the I-mode. 

Figure 2. Stability of a geometrically slim, hot two temperature disc. In (a) no advection 
is included, the solid, long-dashed and short-dashed lines correspond to the acoustic modes, 
thermal mode and viscous mode respectively, (b) shows the case with very little advection 
but without thermal diffusion (q=0.01, f=0). The solid and dotted lines correspond to 
the O-mode and the I-mode. In (c) the differences of thermal instabilities in two kinds of 
optically thin discs are indicated. The heavy lines represent the case of a two-temperature 
disc and light lines represent the case of bremsstrahlung disc. The dot long-dashed lines 
correspond to the thermal modes when little advection is included (q=0.01). (d) shows the 
case similar with (b) but with thermal diffusion included (f=3). The thermal and viscous 
modes are mixed together in this case. 

Figure 3. Stability of a geometrically slim, advection-dominated two temperature disc. 
The solid, dotted, long-dashed and short-dashed lines correspond to the O-mode, I-mode, 
thermal and viscous mode respectively. The parameters (m, f) are (0.01, 0) in (a), (0.1, 
0) in (b) and (0.01, 3) in (c). The thermal and viscous modes are mixed together in (c) 
where the thermal diffusion is included. 
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